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Abstract 

Here, I summarize briefly main results of the studies of density and 
charge fluctuations in multiparticle production in which I was mainly 
involved during last 3 years. 

The solution of QCD equations for generating functions of parton 
multiplicity distributions reveals new features of cumulant moments os- 
cillating as functions of their rank. Experimental data on hadron mul- 
tiplicity distributions in e^e~ ,hh,hA,AA collisions possess the similar 
features. Contrary, the "more regular" models like Xcj>e predict a dif- 
ferent pattern. Evolution of the moments at smaller phase space bins, 
zeros of the truncated generating functions and the singularity of the 
generating function are briefly discussed. Such studies can provide some 
guide to a new representation of multiparticle processes as compared to 
the common Fock representation. 

Also, the large charge fluctuations can be observed if in some events 
pious tend to be produced in coherent or squeezed isospin states. 

1 Introduction and main results 

During the last decade we worked actively on correlations and fluctuations of 
the number of particles produced within small phase space bins inspired by 
the ideas of intermittency and fractality. The factorial and cumulant moments 
were especially helpful in this respect. Here I would like to show that the 
correlations and fluctuations in large reveal new unexpected features in the 
behavior of the moments as functions of their rank. The larger number of 
particles in the total phase space allows to analyze the moments of higher 
ranks compared to small bins. This is somewhat different view of the density 
fluctuations within the available phase space volume. 

Another interesting problem is related to the charge fluctuations in a single 
event. It is closely connected to the so-called Centauro events observed in 
cosmic rays. 

I present from the very beginning both short review of the problem and 
the main results obtained during the last 3 years, leaving their derivation 
and brief discussion for the next section. Those interested in more detailed 
description (and, in particular, in formulae and Figures demonstrating the 
statements and shown at oral presentation) should use the list of references 



for further reading (e.g., the review papers The density fluctuations are 
given by the width (and by higher moments) of the multiphcity distribution, 
while charge fluctuations show e.g. how the share of neutral pions in single 
events declines from its standard value of 1/3. The first problem is treated in 
the framework of QCD, the second one is considered for coherent and squeezed 
isospin states. 

I would like to stress that QCD predicts the distributions of partons 
(quarks and gluons) while in experiment one gets the distributions of final 
hadrons. Therefore no quantitative comparison has been attempted. To do 
that, one must rely on the Monte-Carlo models with some definite hadroniza- 
tion schemes. However, the qualitative features of both distributions are so 
spectacular and remind each other that one is tempted to confirm once again 
that QCD is a powerful tool for predicting new features of hadron distributions 
as well. 

For a long time, the phenomenological approach doniinated in description 
of multiplicity distributions in multiparticle production □. The very first at- 
tempts to apply QCD formalism to the problem failed because in the simplest 
double-logarithmic approximation it predicts an extremely wide shape of the 
distribution in the global phase space □ (i.e. huge fluctuations} that contra- 
dicts to experimental data. Only recently it became possible □ to get exact 
solutions of QCD equations for the generating functions of multiplicity dis- 
tributions which revealed much narrower shapes and aufli a novel feature of 
cumulant moments as their oscillations at higher ranksQa. The similar oscilla- 

tons have been found in experiment for the moments of hadron distributions 
Their pattern differs drastically from those of the popular phenomenological 
distributions El and of the " non^ingular" (however, possessing the asymptotic 
freedom property) A(/)g-model£2l. The QCD inspired Monte-Carlo models de- 
scribe experiment quite wellEil. 

These findings have several important implications El. They show that: 
l.the QCD distribution belongs to the class of non-infinitely- divisible ones. 
Two corollaries of this statement follow immediately: 

a. The Poissonian cluster models (e.g., the multiperipheral cluster model) are 
ruled out by QCD. 

b. The negative binomial distribution (so popular nowadays in phenomenolog- 
ical fits) is not valid. 

2. the new expansion parameter appears in description of multiparticle pro- 
cesses. 

Since this parameter becomes large when large number of particles are 
involved, it asks for the search of some collective effects and of more convenient 
basis than the common particle number (Fock) representation. 



QCD is also successful in qualitative description of evolution of multiplic- 
ity distributions with decreasingjihase space bins which gives rise to notions 
of intermittency and fractality EflEjO. The fluctuations increase in smaller 
bins. However, there are snroja new problems with locations of the minimum 
of cumulants at small bins t£[E3. 

The experimentally defined truncated generating functions possess-arLia- 
triguing pattern of zeros in the complex plane of an auxiliary variable E3EZHi3. 
It recalls the pattern of Lee- Yang zeros of the grand canonicaLpartition func- 
tion in the complex fugacity plane related_.tD phase transitionE3E3 and asks for 
some collective effects to be searched forLilH. At high multiplicities these zeros 
tend to pinch the positive real axis at the singularity position of the gejoerating 
function and their study can reveal the nature of the singularity (seeEllEj) that 
has far-reaching consequences for the theory of multiple production. 

Besides density fluctuations, there could appear the asymmetry between 
neutral and charged pions distributions. The ideas of chiral models, disori- 
ented chiral condensate, coherent and squeezed isospin states could be useful 
in approaching the problem of charge fluctuations. I describe them briefly in 
a separate section. 

2 Some QCD technicalities and results 

Let us define the multiplicity distribution 

oo 

Pn = (^nl^f^n, (1) 

where cr„ is the cross section of n-particle production processes, and the gen- 
erating function 

oo 

G(z) = 5]P„(l + zr. (2) 

The (normalized) factorial and cumulant moments of the P„ distribution are 

_ E„Pnn(n-l)...(n-(7 + l) _ 1 d^G{z) ^ 

{Y.nPnnY ^ {n)'i dz'i 1^=°' ^•'> 



1 d'?lnG(z)| 



K, = ^^-^^l^., (4) 



where (n) — ^„ Pnn is the average multiplicity. They describe full and genuine 
g-particle correlations, correspondingly. 



Here, I consider QCD without quarks i.e. gluodynamics, since quarks do 
not change quahtative conclusions described below. The generating function 
of the gluon multiplicity distribution in the global phase-space volume satisfies 
the equation 

f\xK{x)j^[G{z,Y + In x)G{z,Y + \n{l-x))-G{z,Y)]. (5) 
(^'^ Jo 

Here Y = Iii(p9/Qq),p is the initial momentum, 9 is the angular width of the 
gluon jet considered, p9 = Q where Q is the jet virtuality, Qo =const, 

I'o - (6) 

TT 

as is the running coupling constant, and the kernel of the equation is 

K{x) = ^-{l-x)[2-x{l-x)]. (7) 

The eq.(||) can be solved exactly for fixed coupling constant! and in higher 
order approximations for the running coupling Q. In the last case, the solution 
of this equation in terms of moments looks as 

rr 7g[l-2/^i7 + /^2(g"7'+g70] 

where the anomalous dimension 7 « 70 + 0(7q). The main prediction is the 
minimum of Hq = Kq/Fq at 

^ ^ 0(70) « 5 (9) 



and subsequent oscillations of the ratio Hq at higher q. Let us note that it 
owes to the singular part of the kernel and is absentE3 in more regular theories 
like A(/)i 

While the above results are valid for gluon distributions in gluon jets (and 
pertain to QCD with quarks taken into accountEl), the similar qualitative fea- 
tures characterize the multiplicity distributions of hadrons in high energy re- 
actions initiated by various particles and nuclei. The numerous demonstration 
of it can be found in papersoO. 

Another important feature of the theoretical results is the presence of the 
product 7oq as a new expansion parameter in all the solutions. Formally, 
it vanishes in asymptotics. However, since 70 « 0.48 even at .Z^-peak this 



parameter is large at any rank q and determines the main properties of the 
moments. One can say that the asymptotics is unreachable, in practice. Let 
us remind that this asymptotics is somewhat similar to the negative binomial 
distribution with very wide shape compared to experimental ones (it gives rise 
to Hq = q~'^ while NBD with k = 2 predicts Hq = 2/q{q+ 1)). Since at higher 
ranks one deals with high multiplicity events the product joq indicates that 
for such processes the usual Fock representation is not convenient. Therefore 
one is tempted to look for a more suitable representation for multiparticle 
processes. 

The multiplicity distributions can be measured not only in the total phase 
space (as has been discussed above for very large phase-space volumes) but 
in any part of it. The most interesting problem here is the law governing 
the growth of fluctuations and its possible departure from a purely statistical 
behavior related to the decrease of the average multiplicity in small bins. Such 
a variation has to be connectecl— with the dynamics of the interactions. In 
particular, it has been proposedcJ to look for the power-law behavior of the 
factorial moments for small rapidity intervals Sy 

Fq ^ {5y)-^'^'i'> ((/.(g) >0) {5y^Q), (10) 

inspired by the idea of intcrmittcncy in turbulence. In the case of statistical 
fluctuations with purely Poisson behavior, the intcrmittcncy indices (j){q) are 
identically equal to zero. 

Experimental data on various processes in a wide energy range support 
this idea, and QCD-^provides a good basis for its explanation as a result of 
parton showers E3^Eja. 

Let us turn now to the g-behavior of moments at small bins. The phe- 
nomenon of the oscillations oLcppiulants discussed above reveals itself here as 
well. According to the theoryO^t^l, the first minimum moves to higher ranks at 
higher energies because more massive jets become available. Another corollary 
is that it should shift to smaller values of q for smaller bins at fixed energy 
because the effective value of the anomalous dimension increases due to lower 
effective masses of subjets (since qmin oc 7(7^ )• While former statement finds 
some support in experiment, the second one does not look to be true (probably, 
due to higher order terms in the relation above). Also, one should always keep 
in mind that lower multiplicities in smaller bins prevent from getting higher 
order moments with good enough precision. 

There is another fascinating feature of multiplicity distributions - it hap- 
pens that zeros of the truncated (if the sum in |^ runs up to n = Nmax) 
generating function form a spectacular pattern in the complex plane of the 
variable z. Namely, they seem to lie close to a single circle. At enlarged values 



of Nmax they move closer to the real axis pinching it at some positive value of 
z. 

No QCD interpretation of the fact exists because it is hard to exploit the 
finite cut-off in analytic calculations. The interest to it stems from the analogy 
to the locations of zeros of the grand canonical partition function as described 
by Lee and Yang who related them to possible phase transitions in statistical 
mechanics. In that case, the variable z plays the role of fugacity, and pinching 
of the real axis implies existence of two phases in the system considered. The 
Feynman- Wilson liquid analogy can be used when applying this idea to particle 
production. However, in my opinion, it would be premature to consider this 
phenomenon as a signature of any phase transition in particle collisions. 

In particle physics, it shows up the location of the singularity of the gen- 
erating function. The number of zeros of truncated generating functions in- 
creases and they tend to move to the singularity point when N^ax ~* oo. Since 
it happens to lie close to the origin, it drastically influences the behavior of mo- 
ments (see (H), (H)), and, therefore, determines the shape of the distribution. 
It explains also why the moment analysis is so sensitive to the tiny details of 
this shape. At the same time the stability of the qualitative features of the 
moments for different reactions is very impressive and implies some common 
dynamics. The study of the singularities is at tLe very early stage now (new 
results have been reported at this Workshop inc3), and one can only say that 
the singularity is positioned closer to the origin in nucleus-nucleus collisions 
and it is farthest in e+e^ that appeals to our intuitive guess. 

Let us discuss possible implications of the results obtained. The very 
existence of the new expansion parameter shows that the particle number rep- 
resentation traditionally used in particle physics becomes inadequate for multi- 
particle production. It asks for the search of another approach. The appealing 
example of the analogous situation has been provided by laser physics where 
the coherent state representation is successful. However, in multiple produc- 
tion neither coherent nor squeezed states look quite promising since they do 
not fit above findings about distributions. Probably, their weighted averages 
giving rise to somewhat similar to the negative binomial distribution would 
be suitable. Anyway, the singularity of the generating function can become 
a starting point for further progress in that respect. This is an attractive 
direction for new research. 

To conclude, I would like to stress that, once again, QCD demonstrates 
its power in predicting new features of particle distributions when dealing with 
parton distributions. 



3 Notes on charge fluctuations 



Above, we have discussed the charged multiplicity distributions which deter- 
mine the density fluctuations of charged pions. We assumed implicitly that 
the same is true for neutral pions as well. However, if confirmed, Centauro 
events in cosmic rays imply that the strong charge asymmetry can be observed 
in some (probably, rather rare) events at very high energies. The perturbative 
QCD is unable to solve the problem since charge asymmetry appears at the 
hadronization stage only. There is no strong charge asymmetry in the com- 
mon models. Therefore, one should rely on different dynamics. The analogy 
to photon states can be useful. 

It is well known that the soft photons are produced in a coherent state. 
If the soft pions tend to be produced in a coherent (or squeezed due to the 
non-linfiax interactinn) state and respect the isospin conservation, then one can 
showEJc3 (see alsoB) that the strong charge asymmetry should be seen in the 
individual events with a noticeable fraction of them having no charged pions at 
all. This is done by considering the projections of these states on the states with 
the definite value of the isospin. Due to isospin conservation, the total isospin 
of the pion system created at high energy is strongly restricted and should be 
much less than the number of pions. The differential distribution of the ratio 
/ of the number of neutral pions to the total number of pions decreases very 
slowly as f~^^^- The similar situation is typical for the chiral model used for 
descjribing the disoriented chiral condensate and the "Baked Alaska" scenaria 
^O. These approaches favor the creation of squeezed isospin states as well 

However, neither coherent nor squeezed states are able to reproduce the 
pattern of moments oscillations described above. No more room is available 
here to discuss thcj)roblem. For more details about coherent isospin states, I 
refer to the papered presented at this Workshop. 

I should apologize that in these notes I described briefly the main state- 
ments about density and charge fluctuations in particle production leaving 
aside numerous details contained in the papers presented in the list of refer- 
ences. However, the size of the presentation prevents the detailed review. 

Acknowledgment s 

I am grateful to W. Kittcl for inviting me to participate in this workshop 
and for financial support. 

This work is partly supported by the Russian Fund for Fundamental Re- 
search (grant 96-02-16347a) and by INTAS grant 93-79. 



References 



1. I.M. Dremin, Physics-Uspekhi, 37 (1994) 715. 

2. I.M. Dremin, and A.V. Leonidov, Uspekhi Fiz. Nauk, 165 (1995) 723; 
Physics-Uspekhi, July 1995. 

3. E.A. De Wolf, I.M. Dremin, and W. Kittel, Phys. Rep. 270 (1996) 1. 

4. Z. Koba, H.B. Nielsen, and P. Olesen, Nucl. Phys. B40 (1972) 317. 

5. Yu.L. Dokshitzer, V.A. Khoze, A.H. Mueller, and S.I. Troyan, Basics 
of Perturbative QCD, Gif-sur-Yvette, Editions Prontieres, ed. J. Tran 
Thanh Van, 1991. 

6. I.M. Dremin, Phys. Lett. B313 (1993) 209. 

7. I.M. Dremin, and V.A. Nechitailo, JETP Lett. 58 (1993) 945. 

8. I.M. Dremin, and R.C. Hwa, Phys. Rev D49 (1994) 5805. 

9. I.M. Dremin, V. Arena, G. Boca et al, Phys. Lett. B336 (1994) 119. 

10. I.M. Dremin, and V.A. Nechitailo, Yad. Fiz. (Phys. Atom. Nucl.) (to 
be published). 

11. A. Capella, I.M. Dremin, V.A. Nechitailo, and J. Tran Thanh Van, Z. 
Phys. (1996) (to be published). 

12. Yu.L. Dokshitzer, and I.M. Dremin, Nucl. Phys. B402 (1993) 139. 

13. W. Ochs, and J. Wosiek, Phys. Lett. B305 (1993) 144. 

14. Ph. Brax, J.L. Meunier, and R. Peschansky, Z. Phys. C62 (1994) 649. 

15. I.M. Dremin, Phys. Lett. B341 (1994) 95. 

16. I.M. Dremin, in Proc. 4th Int. Workshop on Relativistic Aspects of 
Nuclear Physics, Rio de Janeiro, Aug. 1995, WSPC (to be published). 

17. E.A. Do Wolf, A note on multiplicity generating functions in the complex 
plane (unpublished); Proc. of 24th Int. Multiparticle Dynamics Sympo- 
sium, Vietri, 1994, WSPC, Eds. A. Giovanini, S. Lupia, A. Ugoccioni, 
p. 15. 

18. G. Gianini, in Proc. of 25tli Int. Multiparticle Dynamics Symposium, 
Slovakia, Sept. 1995, WSPC, Singapore (to be published). 

19. C.N. Yang, and T.D. Lee, Phys. Rev. 87 (1952) 404. 

20. T.D. Lcc, and C.N. Yang, Phys. Rev. 87 (1952) 410. 

21. I.M. Dremin, JETP Lett. 60 (1994) 757. 

22. G. Gianini, the talk at this Workshop. 

23. A. Bialas, and R. Peschanski, Nucl. Phys. B273 (1986) 703. 

24. I.V. Andrccv, JETP Lett. 33 (1981) 367. 

25. I.M. Dremin, and R.C. Hwa, Phys. Rev. D53 (1996) 1216. 

26. A.A. Anselm, Phys. Lett. B217 (1989) 169. 

27. J.D. Bjorken, Int. J. Mod. Phys. A7 (1992) 4819. 

28. I.I. Kogan, Phys. Rev. D48 (1993) 3971. 

29. I.V. Andreev, the talk at this Workshop. 



